Quantum integrability and nonintegrability in the spin-boson model 
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We study the spectral properties of a spin-boson Hamiltonian that depends on two continuous 
parameters < A < oo (interaction strength) and < a < n/2 (integrability switch). In the classical 
limit this system has two distinct integrable regimes, a — and a — n/2. For each integrable regime 
we can express the quantum Hamiltonian as a function of two action operators. Their eigenvalues 
(multiples of h) are the natural quantum numbers for the complete level spectrum. This functional 
dependence cannot be extended into the nonintegrable regime (0 < a < n/2). Here level crossings 
are prohibited and the level spectrum is naturally described by a single (energy sorting) quantum 
number. In consequence, the tracking of individual eigenstates along closed paths through both 
regimes leads to conflicting assignments of quantum numbers. This effect is a useful and reliable 
indicator of quantum chaos - a diagnostic tool that is independent of any level-statistical analysis. 

PACS numbers: 05.45.-a; 05.45.Mt 



I. INTRODUCTION 

Classical integrability of a system with two degrees 
of freedom guarantees that the Hamiltonian can be ex- 
pressed as a piecewise smooth function of two action 
coordinates: H(pi,qi]p2,(h) = Hc{Ji,J2)- No such 
functional relation exists if the system is nonintegrable 
[2 0, B 0] ■ Geometrically speaking for a parametric sys- 
tem with parameters subject to an integrability condi- 
tion, there exist complete foliations of invariant tori in 
phase space for all parameter points in the integrable 
regime. Throughout the nonintegrable regime the foli- 
ation is partially destroyed. Some tori are replaced by 
chaotic trajectories, cantori, and unstable periodic tra- 
jectories. The surviving tori in the nonintegrable regime 
are no longer dense anywhere in phase space, but each 
one of them can still be characterized by two local action 
coordinates Ji, J 2 via line integrals § Pidqi along pairs of 
topologically independent closed paths. 

In the integrable regime, the natural label of an invari- 
ant torus is its set of action coordinates ( J±, J2). Tracking 
an invariant torus along a path through the integrable 
regime of parameter space means that we observe how 
the torus with this specification changes its position and 
shape in phase space. In the nonintegrable regime, where 
all intact tori are separated by chaotic phase flow, an in- 
dividual torus can no longer be characterized by fixed 
values of J\ , J2 ■ Tracking a surviving invariant torus 
along a path through the nonintegrable regime of param- 
eter space now means that we observe it in isolation from 
other tori. The local action coordinates vary smoothly as 
the torus changes its location and shape in phase space. 

Now let us attempt to track one torus along a closed 
path that lies partly in the integrable regime and partly 
in the nonintegrable regime, assuming that it does sur- 
vive the presence of chaos. Inside the integrable regime 
the identity of the torus is determined by constant val- 
ues of the action coordinates, while outside that regime 
the action coordinates vary with the shape of the iso- 



lated torus. The values of J\, J2 at the end of the closed 
path will, in general, be different from the starting val- 
ues, implying that the individuality of a torus cannot be 
maintained. No loss of individuality is suffered by tori 
along closed paths embedded in the integrable regime or 
for surviving tori along closed paths in the nonintegrable 
regime [5|. 

There exists a quantum counterpart to this 'crisis of 
identity' as will be demonstrated. It can be employed to 
discriminate between regimes of integrability and nonin- 
tegrability on purely quantum mechanical grounds. Here 
we show the workings of this diagnostic tool in the con- 
text of the spin-boson model @, 0, H, Q , 

H = fvujBd 1 a + h^sSz 
+ A cos a (S + a + S-a^) 
+ A sin a (S+a^ + S-a) , (1) 

one of the simplest nontrivial models describing nonrel- 
ativistically the interaction between an atom and a ra- 
diation field This model has also been used to de- 
scribe the interaction between electronic and vibrational 
degrees of freedom in molecules and solids. The relation 
between classical and quantum integrability of (ji} has 
been the object of previous investigations [a, fllj. 

The first two terms in (fTJ) describe one mode of the 
electromagnetic field and a (2a + l)-lcvcl atom, respec- 
tively. The coupling between the two degrees of freedom 
has strength A and depends on a continuous parameter 
a that connects two regimes for which this model is in- 
tegrable in the classical limit. The classical integrability 
for a = and a = n/2 is established by a second integral 
of the motion. The case a = is known as the rotating 
wave approximation in quantum optics. Early studies in 
one or the other classical limit of the spin-boson model 
revealed chaotic phase space flow turning regular in the 
rotating wave approximation 

In the two-dimensional parameter space spanned by 
the (polar) coordinates (A, a), the two integrable regimes 
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are located on two perpendicular straight lines that in- 
tersect each other at the point of zero coupling strength. 
Each quadrant of this parameter plane represents a non- 
integrable regime. Henceforth we consider the parameter 
range 0<A<oo, 0<a< n/2. 

In preparation of our main theme we first discuss the 
classical integrability condition of the spin-boson model 
(Sec. HIl and then the classification of its quantum energy 
levels fSec. HIip and certain quantum invariants fSec. IIV[) 
by distinct sets of quantum numbers in the integrable 
and nonintegrable regimes. This distinction has a deeper 
meaning, which we will further discuss in Sec. |Vl and 
which we will employ in Sec. IVII for the identification of 
the two regimes in purely quantum mechanical terms. 



II. INTEGRABILITY CONDITION 



In taking the classical limit h — > 0, a — * oo of the 
spin-boson model, we renormalize the coupling constant, 
A = (ft/2) 3 / 2 A, substitute 



Mujb 
2h 



1 



2hMuj B 



ip 



(2) 



for the boson operators and, via hy a{a + 1) = s, convert 
the spin-CT operator into a classical 3-component vector 
of fixed length: 



(S x , S y , S z ) = s (sin i? cos tp, sin d sin tp, cost)). (3) 



dp/dt 



-dH/dx, and dS/dt = -S x dH/dS read 
p A 



M 2^/MU^ 



5,, (sin a — cos a), 



2 AV Muj b „ , 
p = —Mu) B x 6 x (cosa + sin a), 



S x — —u>sS. 
S y = ujsS x 



Ap 



" 2^/A1Ub~ 
Ax\/ MuJ B 



S z (cos a — sin a) , 
S z (cosa + sin a), 



S z = 



S y (cos a + sin a) 



Ap 



2y/Mu7^ 



S x (cos a — sin a) 



(5a) 

(5b) 
(5c) 

(5d) 
(5e) 



The phase flow generated by these equations is, in gen- 
eral, chaotic. Chaos gives way to a fully intact torus 
structure at a = 0, ir/2. The integrability of these cases 
is established by the fact that one or the other of the two 
functions. 



I = 
K = 



1 



2Mu B ' 2 MUJBX2 + Sz > 



Mu B x 2 -S z , 



2Muj b 2 
whose time evolution is determined by [l4| 



/ = {HJ} = Asina ( f Sx - JMuj b xSA 
K = {H,K} = A cos a \ -fj= + \/Mu> B xS v 



(6a) 
(6b) 

(7a) 
(7b) 



becomes a second integral of the motion. The case a = 
is known as the Jaynes-Cummings model [l5T |. The im- 
pact of the classical integrability conditions on the quan- 
tum system is the main theme of this study. 



The spin-boson Hamiltonian (fTJ) thus turns into the en- 
ergy function of two linear one-degree-of-freedom systems 
- a harmonic oscillator and a classical spin in a constant 
magnetic field - with a nonlinear coupling: 



H = ^ + -Mulx 2 



2M 
+ - A cos a 



-A sin i 



■u) S S z 



\J Mu B xS x 



1 



\J Muj b xS x 



yj Mui B ^ V 
1 



y/Mu)B 



pS, 

pSy 



(4) 



A set of canonical coordinates is (p, x; scosf?, tp). 
The equations of motion for the physical variables 
(x,p, S x , S y , S z ) inferred from (f5|) via dx/dt = dH/dp, 



III. ENERGY LEVELS 

For the analytic or numerical solution of the spin- 
boson model (H|), it is convenient to use the prod- 
uct vectors of the noninteracting system, \m,n), m 

0. 1.2 2(7, n = 0, 1, 2, . . ., as a basis. The relevant 

operators act on this basis as follows: 



(<t - S z )\m,n 
S+\m, n 
S- 1 m, n 
a' |m, n 
a|m, n 



= m 1 77i, n), (8a) 
= y / m(2(T - m + l)\m - 1, n), (8b) 



y/(2a - m)(m + l)\m + 1, n), (8c) 

(8d) 
(8e) 



y/n + l\m, 77 + 1) 



The Hamiltonian matrix can thus be assembled from the 
diagonal elements 

(m, n\S z \m, n) — a — m, (m,n\a' a\m,n) = n, (9) 
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n 



1 










2 
m 



2a 



FIG. 1: Basis vectors \m, n) with positive parity (full squares) 
and negative parity (open squares) as coupled by matrix el- 
ements of operators S+a,S-a* (solid lines) and S+a^,S-a 
(dashed lines) of Hamiltonian |T} 



and from the off-diagonal elements 



(m, n\S+a\m + 1, n + 1) = \J (2a — m)(m+ l)(n + 1), 



(m, n\S-a)\m — 1, n — 1) = \/ (2er + 1 — m) 



mn, 



(m, n|5'+a^|m + l,n — 1) = \/ (2<r — m)(m + l)n, 



(m, n\S_a\m — l,n + 1) = y (2er + 1 — m)m(n + 1). 

The structure of this matrix is illustrated in Fig. [T] The 
solid lines represent matrix elements generated by the 
first interaction term in ([1]), and the dashed lines repre- 
sent matrix elements which arise in the second interac- 
tion term. The two-sublattice structure is a reflection of 
parity conservation. The parity operator, 



P=(-l)° 



(10) 



commutes with H for arbitrary A, a. It divides the eigen- 
states into two symmetry classes. States with P = +1 
(P = —1) involve basis vectors with even m + n (odd 
m + n) only. 

If a — only the solid bonds are present and if a = n/2 
only the dashed bonds. In either case the Hamiltonian 
matrix is reduced to invariant blocks of size 2<r + 1. If 
< a < 7r/2 wc must deal with an infinite matrix. In this 
study we set cos = = w except where indicated other- 
wise. In the following we analyze the level spectrum for 
various cases at a — 0, ir/2: systems with a = |, a = 1 
for arbitrary n, and systems with arbitrary a asymptot- 
ically for large n. 



A. Spin a = | 

The integrable spin-boson model with a = h involves 
only 2x2 matrices. If a — 0, the eigenvectors happen to 
be independent of the interaction strength, 



l^i.o) = |1,0) 



(11a) 



|Vi,n> = -7= {|1,«> + 10,71-1)}, n=l,2,..., (lib) 
l^o.n) = -j= {|l,n + 1) - |0,n)} , n = 0, 1, . . . , (11c) 

and the energy eigenvalues (for n = 0, 1, 2, . . .) are 

E 1>n = Hu)(n-l 12) + Ay/n, (12a) 
E , n = huj(n + 1/2) - AVn + 1. (12b) 
If a = 7r/2 the eigenvectors do depend on A: 

|Vi,o) = |0,0) (13a) 
iV'i.n) = ao,n\0, n) +bo, n \l,Ti— 1), n=l,2, (13b) 
iV'o.n) = ai,7i|0,n+ 1) - &i, n +i|l,n), n = 0,1,... (13c) 
with 



«o, 



%/A~T 



bo.n — 



fx 



0£) 



^(Ai-VA;)' 



O-l, n-l 



01,n-l — 



a/2(Ai + yfK) ' 

a/AT+1 



\/2(Ai + V%)' 



where A = (A/hco) 2 , A„ = 1 + nA. The associated energy 
eigenvalues (for ra = 0, 1, 2, . . .) are 

Ei tn = hco(n - 1/2) + bi/A^, (14a) 
£ ,n = + 1/2) - tuo y/X n+1 . (14b) 

B. Spin a = 1 

The case a = 1 at integrability involves the solution of 
cubic equations. Here we list the (A-independent) eigen- 
vectors and the associated energy eigenvalues for a — 0. 
We have |^, ) = |0,0), 1^,0 = (|1,0) + |0, l))/y/2, 
|^2,i) = (|1,0) - |0,1))/V2, with energies E lfi = -hu>, 
E\^i = v2A, E2.1 = — v2A, respectively, and for n > 2 
the results are 



IV>2,n) 

\i>a,n) 



\n- 


1 


An - 


2 


' n 


An - 


2 


1 n 


2n- 


1 


1 n — 


1 


An - 


2 


1 n 



An- 2 



|0,n-2) 

|2,n), 

\0,n-2) 

\0,n-2) 

|2,n), 



V2 



l,n-l 



n- 1 
2n - 1 

1 



12, n) 



l,n-l 



(15a) 
(15b) 

(15c) 
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with energies 



Ei n = fku(n — 1) 
E 2 , n = ^>{n - 1), 
£3,71 = hw{n - 1) 



Ay/4r 



AVTr 



(16a) 
(16b) 
(16c) 



C. Arbitrary Spin a 

A simple analytic solution exists for arbitrary a in the 
asymptotic regime of large n. Consider the (2a + 1)- 
dimensional invariant block of (fTJ) at a = formed by 
the basis vectors \2a — to, n — to), to = 0, 1, . . . , 2a. It is 
tridiagonal with elements 

(2a — m,n — m\H\2a — m,n — to) = fku(n — a), 
(2a — to, n — m\H\2a — to, n — to — 1) 



k\/2a(n - m). 



For n 3> a we can write 



H = hiu(n - a)£ + 2ky/nS x + O [ -= , 



(17) 



where £ is the (2a + l)-dimensional unit matrix and S x is 
the irreducible representation of the spin operator with 
the same dimensionality. The asymptotic eigenvalues of 
this matrix are 



E mM ~ huj(n — a) + 2h^fn(a — to) 



(18) 



for to = 0, . . . , 2a. The corresponding analysis carried 
out for a — tt/2 yields the matrix 



H = hw(n - a)£ + 2hwS z + 2Ay^S x + O 



(19) 



with asymptotic energy eigenvalues (for to = 0, ... , 2a) 



hw(n — a) + 2hio\] X n (a — to) 



(20) 



which indeed commute with (fTJ under exactly the same 
conditions as in the classical limit. We have 

[H,I] = 2Asma(S_a- 5 + a f ), (22a) 
[H, K] = 2Acosa(5+a-5_a t ). (22b) 

However, quantum integrability cannot be inferred from 
quantum invariants as simply as classical integrability 
can be inferred from integrals of the motion (analytic 
invariants). Commuting operators can always be con- 
structed irrespective of whether the model is (classically) 
integrable or not [l6|, H3] ■ The parity operator (fTT)|) , for 
example, which can be expressed as a function of cither 
invariant I 01 K, 



p _ e i7r(//fi— er) _ e i-n(K/h+o) 



(23) 



commutes with H for arbitrary a. More generally, any 
operator A that is not already an invariant, [H, A] =/= 0, 
can be turned into an invariant via time average. In the 
energy representation, the time average strips A of all its 
off-diagonal elements. The resulting operator 1a = (A) 
thus commutes with H by construction [3 [l!| . 

The fact is that in the classical limit neither the parity 
operator nor any of the artificially constructed quantum 
invariants will turn into analytic invariants (integrals of 
the motion) if the phase flow is chaotic. Such quantum 
invariants either lose their meaning altogether or turn 
into nonanalytic invariants (l9l [20| . 

The distinctive attributes of quantum invariants in the 
integrable and nonintegrable regimes of a quantum sys- 
tem are subtle but not ambiguous. Here we use 



I A = (A), A = a^(S-+S+). 



(24) 



For a = i, its eigenvalues at a = can be calculated 
from the eigenvectors (jlip , 



{Ah 



1 



^ (A) , 



1 



Vn+T, 



(25) 



2 V ' x ' " 2 
and its eigenvalues at a = tt/2 from the eigenvectors (fT3 



Note that in all cases pertaining to the integrable 
regimes a = 0ora = 7r/2 the energy levels are naturally 
labelled by the two quantum numbers to, n. The par- 
ity becomes P = (— l)' Tl+ ™. In the nonintegrable regime 
< a < tt/2, by contrast, the numerical analysis suggests 
the use of a single (energy sorting) quantum number k 
for all levels of given parity. 



IV. QUANTUM INVARIANTS 

The quantum counterparts of the two analytic invari- 
ants ([6|) are the operators 



(A)l.n 
(A)o, n 



(A n 

2(X n ~ %/An) 

(K1+1 — K1+1 



1) 



2(A n + l + y/ A n +l) 



(26a) 
(26b) 



For a = 1 and a = 0, we obtain the following eigenvalues 
from the eigenvectors (|15l) : 



{Ah 



-(A) 



3,71 



y/n-1/2, (A) %n = 0, (27) 



Asymptotically for large n, we can evaluate the eigenval- 
ues for arbitrary a. The results for a = read 



H(a^a + S z ), K = h(a)a- S z ), 



(21) 



(A)r 



(a — to) y/n, 



(28) 
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and for a — tt/2 we have 



(a — m)i/n 
y/T+tfnX' 



(29) 



Numerical results of (A)/, for 0<a<7r/2of systems 
with a = i, 1, | were reported previously [9j. The pat- 
terns of points (E m>n , (A} mtTl ) for integrable cases were 
found to be strikingly different from the pattern of points 
(Ek, (A)k) for nonintegrable cases. Here this difference 
will be used as a demarcation tool for regimes of integra- 
bility and nonintegrability. 



Quantum mechanically, the interaction distorts the 
eigenvalue spectrum and modifies the selection rules of 
transition rates. Quantum properties that are as sensi- 
tive to the integrability status as their classical counter- 
parts do exist and have previously been ex plor ed in the 
context of a different model system [2l], I22I |23|| . 

These properties are directly related to the existence 
of action operators as constituent elements of the Hamil- 
tonian such as discussed in Sec. IV Al for the noninter- 
acting system. In the interacting cases, the existence of 
action operators can again be demonstrated directly for 
a = 0, 7r/2, and their nonexistence for < a < tt/2 can 
be demonstrated indirectly. 



V. QUANTUM ACTIONS 

One hallmark of integrability in a quantum system 
with two degrees of freedom is that the Hamiltonian can 
be expressed as a function of two action operators J\ , J2 , 
i.e. of two quantum invariants whose spectra consist of 
equidistant levels. 



A. From A = to A > 

In the absence of the spin-boson interaction (A = 0), 
the two action operators are 



Ji = h(a — S z ), Ji = Ha* a 

with integer eigenvalues (in units of h) 

Ji = mh, m — 0, 1, ... , 2a, 
J-2 = nh, n = 0, 1, . . . , 



(30) 



(31a) 
(31b) 



as in ^j. The Hamiltonian, Ho = hujBCva + hu>sS z , and 
the two quantum invariants (|21[) are expressible as linear 
combinations of J\ , J2 . 

Classically, the contribution of each degree of freedom 
to H = p 2 /2M + \Mlo\x 1 +lo s S z is transformed into a 
function of one action coordinate by a separate canonical 
transformation: (S z , if) — » (Ji, 8\) with S z = s— Ji, (p = 
—81; and (p, x) — > (J2,6*2) with p — y/2J 2 MujB cos 82, 
x = y / 2J2/AluiBsm82- The transformed Hamiltonian 
and the two classical invariants ([6]) are linear functions 
of Ji, J 2 just as in quantum mechanics. The exact quan- 
tum spectra of H ,I,K can then be recovered exactly 
via semiclassical quantization, i.e. by substituting the ac- 
tions quantized according to (|3 1 [) into the classical Hamil- 
tonian. 

Classically, the interaction renders the equations of 
motion, Eqs. ([5]), nonlinear. However, the effects of an- 
harmonicity in the time evolution depend sensitively on 
whether integrability is sustained or destroyed by the in- 
teraction. Integrability for a — 0,7r/2 dictates that the 
phase flow is exclusively toroidal. For < a < tt/2 
chaotic phase flow is omnipresent albeit constrained by 
surviving tori. 



B. a = |, a = 



The unitary transformation which diagonalizes the 
Hamiltonian ([T]) for a = i and a = 0, expressed in terms 
of spin and boson operators, reads 



U A = Pq + -i ( -25, + -^a^S- + aS +n L= 
V2 \ va'a Va*a 



(32) 

where P Q A = |1,0)(1,0|, = 1- |0,0)(0,0| - |1,0)(1,0|. 
The operators 



T z = V A S Z V A X - P^S Z - ~G? 

tfb = Ua^clUa 1 = a f a - S Z P A + ^G A 



with 



G A = aS^ 



(33a) 
(33b) 

(34) 



are diagonal in the energy representation: 

T z \ipm,n) = (f - m)\ip m ,n), (35a) 
fr^lVw) = n\ip m , n )- (35b) 
Hence the quantum actions with eigenvalues pip are 

Ji = h{a - T z ), J 2 = htfb. (36) 
Applying Ua to the Hamiltonian yields 
UaHUa 1 = f^u{b^b + T z ) 
+ A 



(37) 



which, together with (|3"6"|) , describes the functional rela- 
tion between H and J\ , J 2 ■ 



C. a= \, a = tt/2 

The same method also produces the quantum actions 
for the integrable case s = |, a = tt/2 of the spin-boson 
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Hamiltonian ([1} . Here the block-diagonal unitary trans- 
formation, Ub to be used can also be expressed in terms 
of spin and boson operators but has a more complicated 
structure than U a- The operators 

T z = UbSzUb 1 = ip<f + (G B ))Q*, (38a) 

6+6 = U B <JaU B = (a f a- S z + G B ) Qo , (38b) 

with G B = Gf + Gf + Gf, 

of- 1 



^16(l + Aat a ) ^/16 (1 + A (aJa + 1)) ' 



G 



Go 



5 Z 



5'. 



v/4 (1 + Aoto) v' 4 ( 1 + A ( ata + 1 ))' 
1 + 2S 2 f 1 + 2S 2 



V64A + at a ' 



4V64A + at a 

and P<f = |0,0>(0,0|, Q B = 1 - |0,0)(0,0|, again satisfy 
(f3"5)) and are related to quantum actions via (|3U|) . The 
functional dependence of the transformed Hamiltonian 
on the actions is different from (l37l): 



1 -L 97 1 

UsHUg 1 = hu(tfb - T z ) + — + A6tfo 



1-2T 2 



v /l + A(6t6+l). (39) 



D. 



cr > 



The results of Sees. IV Bl and IV CI are generalizable to 
arbitrary tr, albeit for the price of a higher and higher 
calculational effort. The case a = l,a = can still be 
presented compactly. The unitary transformation Uc to 
be used in this case is now determined by the eigenvectors 
(fT5|) and yields 



T z = U C S Z U C Y = 
6+6 = Uc^aUc 1 
= (ata + S z -\ 



Gf) gf 



where 



Gf 



S z Sz 
2\/4ata - 2 
1-S 2 



a+S- 



Gf = 1 



V4a+a + 2 
1-52 



a+S- + 
05*+ + 



2^2 Pi , 



l-5 7 2 
-I — z 

V4ata + 2 
S z + 5* 2 
V4ata + 6 



(40a) 
(40b) 



aS- 



+ 



\/2 2^ 

• c = |0,0)<0,0|,Pf = |1,0>(1,0| + |0,1)(0,1|,QC = 
c _ pC rpj :ic transformed Hamiltonian becomes 



1 — P° — P l 

1 M) r \ 



UcHUc 1 = hu (6+6 + T 2 



A 
71 



[(3T 2 2 - - 2) Pj 



c 



+ ( (T 2 - T z )y / 26t6-l- (T 2 2 + T z )y/2tfb + 3) Qf 



(41) 



Pa and Pc are special cases for a = | and a = 1, 
respectively, of a unitary transformation P 1 (cr, A) that 
diagonalizes |T]) at a — for arbitrary values of a. This 
transformation turns out to be A-indcpcndcnt for the two 
cases we have worked out. It may well be A- independent 
for arbitrary a. Likewise, Ub is the special case for = \ 
of a unitary transformation P 2 (cr, A) that diagonalizes 
(fT|) at a = 7r/2 for arbitrary a. That transformation is 
manifestly A-dependent. 

The end-product of these unitary transformations are 
two functions Hq\t z , 6+6; A) = H Q X) {Jx, J 2 ; A) and 

P^ 2) (T 2 ,6t6;A) = H Q 2) (J 1 ,J 2 ;A), which express the 
functional dependence of the Hamiltonian on action op- 
erators in the two integrable regimes a = and a = tt/2, 
respectively. The leading term of an asymptotic expan- 
sion at high boson occupancy and unrestricted spin state 
of these functions can be inferred from (TIT)) and (|T5|) : 

H Q = huj(tfb±T z ) + Av / 6+6 + 0(l), (42) 

where the operators P 2 ,6+6 again satisfy (j35|) and the 
upper (lower) sign pertains to a = (a = tt/2). 

We expect a semiclassical regime to exist at large 
spin and/or boson quantum numbers where the functions 

(2), 



q\Ji, J2', A) and Pq ; (P,J 2 ;A) connect with func- 



tions Hq(Ji, J2', A) and Hq(J\, J 2 ; A) of classical ac- 



tions. However, the identification of the semiclassical 
regime requires a complete solution of the classical equa- 
tions of motion (O, a task still outstanding. 

The connections between the quantum and classical 
functional dependences of Hamiltonian on actions was 
investigated in a previous study for an integrable two- 
spin model and for the (integrable) circular billiard model 
[231 ] . There we found subtle quantum effects that restrict 
the range of the semiclassical regime in unexpected ways. 
That may also be the case in the spin-boson model. How- 
ever, the point we wish to emphasize in this study is a 
different one. 



VI. TRACKING EIGENSTATES 

The goal is to demonstrate that the functions 
Hq\J\, J 2 ; A) and Hq\j±, Ji\ A) cannot be extended 
in any consistent way into the region of nonintegrability 
in the (A, a)-plane. The functions Hq^ and Hq make 
it possible to label all eigenstates of {I} by the two ac- 
tion quantum numbers ra,n as defined in (|3"Tj) and to 
track them with no ambiguity through each one of the 
two integrable regimes a = and a — tt/2. The non- 
extendability of the two functions Hq and Hq into a 
function Hq(J\, J 2 ; A, a) translates into the impossibility 
of consistently assigning action quantum numbers m, n to 
the eigenstates in the entire parameter range < A < 00, 
< a < tt/2. 

One way of keeping track of eigenstates \v) of (fT|) is 
to determine how the eigenvalues of quantum invariants 
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vary along some path in the (A, a)-plane. For the pur- 
pose of this demonstration, we focus on the eigenvalues 
(H) v = E v of the Hamiltonian {]]) with spin quantum 
number a = i and the eigenvalues (A) v of the quantum 
invariant I a as defined in (|24[) . 

Within each of the two integrable regimes, both sets 
of eigenvalues have an explicitly known (discrete) depen- 
dence on the action quantum numbers v = (m, n) and 
an explicitly known (continuous) dependence on the in- 
teraction strength A. The functional relations are stated 
in Eqs. (HJ), (jUJ) for a = and (14]), ([26]) for a = tt/2. 



A. Level Crossings 

In panels (a) and (b) of Fig. [2] we have plotted one 
quantum invariant versus the other for all states with 
positive parity up to a certain energy. In both panels we 
observe two vertically displaced rows of states. States in 
the top and bottom rows have action quantum numbers 
(l,n) and (0,n), respectively. 

The observed arrangement of states is due to the fact 
that {A) m ,n ~ Vn but E m ^ n ~ n in leading order. No- 
tice that the spacings between successive energy levels in 
each row vary slowly, and at different rates in the top 
and bottom rows. To enhance the visibility of this effect 
we have connected successive energy levels in each panel 
by dashed lines. The spacings are somewhat larger in the 
top row compared to the bottom row, causing instances 
in both panels where two consecutive states of the bottom 
row fit into the space between two states of the top row. 
These instances where the alternating (top/bottom) se- 
quence is broken mark locations where energy levels from 
opposite rows can fall arbitrarily close to each other. 

When we increase the interaction strength A gradually, 
the states in the top row of Fig. [Ua) move toward the 
right and the states in the bottom row toward the left. 
The same observation can be made in Fig. H^b). Here 
the shift also contains a small vertical component. We 
have singled out one pair of nearly degenerate states in 
Fig. [Ha) an d another pair in Fig. [2(b). Each pair is 
marked by full circles. In panels (a) and (b) of Fig. [3] 
we have plotted the traces of these states in the plane 
of invariants as the interaction strength is increased by a 
certain amount. 

The gradual change of A causes a cascade of level cross- 
ings between states from opposite rows. For the two pairs 
of tagged states, the crossings occur at the point marked 
by an asterisk on their traces. States from opposite rows 
undergo level crossings even though they have the same 
parity. What matters are the functional relations Hq^ 

(2) 

and Hq ' established previously. They remove any possi- 
ble cause for level collisions (avoided crossings) between 
states from opposite rows as they move (energetically) in 
opposite directions when A is increased. 



B. Level Collisons 

A very different scenario unfurls when we plot the two 
quantum invariants for a nonintegrable case. What hap- 
pens when we change the integrability parameter from 
a = [Fig. Ufa)] or from a = tt/2 [Fig. [^b)] to a = tt/4 
is illustrated in Fig. [5(c). Here the states that used to live 
in different worlds (top row with action quantum num- 
ber m = 1 and bottom row with m = 0) now suddenly 
get into each other's way. Since they are prohibited from 
undergoing any level crossings, it is now appropriate to 
label them by the energy sorting quantum number k. 

In those parts of the spectrum where the energy level 
spacings are large, the loss of integrability has no visible 
effect on the quantum invariants. That is the case near 
the left and right border areas of Fig. [5(c). Here the 
two rows of states remain largely intact. However, near 
the center of the panel, where small energy level spac- 
ings occur, the eigenvectors of nearly degenerate levels 
affect each other strongly. The most conspicuous effect 
is a strong vertical displacement of the two states from 
the row positions toward each other. Less conspicuous in 
Fig. [3](c) but of even greater importance is the small hor- 
izontal displacement of the two nearly degenerate states 
away from each other. The effect of nonintegrability is 
that energy levels exert a short-distance repulsion on each 
other. At the same time, expectation values in general 
and the quantum invariant (A) j. in particular tend to be- 
come less differentiated than they were in the integrable 
case. 

When we again increase the interaction strength A, 
now at fixed a = 7r/4 in the nonintegrable regime, we find 
that no levels with equal parity ever undergo a crossing. 
As in the integrable cases, the states with (A)k > have 
a tendency to move toward the right and the states with 
(A)k < toward the left. 

Inevitably, these trends put states on opposite sides 
of {A)k = on a collision course. When two such states 
approach one another, the state starting out with (A)k > 
swings down as it moves to the right and the state 
with (A)k < swings up as it moves to the left. The 
two states reach their closest energetic approach when 
their vertical positions are about the same. After that, 
the state coming from below continues its upswing, but 
now it is moving to the right to join the right-moving 
upper row of states. Meanwhile, the state coming from 
above continues its downswing to join the left-moving 
lower row of states. One such level collision, between the 
tagged states in Fig. (He), is shown in Fig. [3(c). 

C. Quantum Numbers in Conflict 

Looking at the spectrum of the spin-boson model ([1]) 
in the plane of invariants {E V ,{A) V ) as the interaction 
strength A increases gradually, reveals strikingly different 
patterns of coordinated motion of all states with given 
parity, depending on whether the parameter a is set to 
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FIG. 2: Quantum invariant {A) v = (u\a\S- + S+)\u) versus 
quantum invariant E v — (u\H\u) over some energy range for 
the eigenstates \v) with parity P = +1 of the spin-boson 
model Q with a = |, fkd = 1, X = (A/hco) 2 = 0.09, and 
(a) a = 0, (b) a = 7f/2, (c) a: = 7r/4. In the integrable 
regimes we use the action quantum numbers v = (m, n) and in 
the nonintegrable regime we use the energy sorting quantum 
number v — k. One pair of states in each panel (full circles) 
is tagged for further use in Fig. [3] 



an integrable regime (a — 0, 7r/2) or fixed within the 
nonintegrable regime (0 < a < 7r/2). 



FIG. 3: Trace of one pair of eigenstates \v) with parity 
P = +1 (identified by full circles in Fig. [2} in the plane of 
quantum invariants (E v , {A) v as the interaction parameter A 
is increased a specified amount at constant value (a) a = 0, 
(b) q = 7r/2, (c) q = 7r/4 of the integrability parameter. In 
the integrable regimes we use v = (m, n) and in the noninte- 
grable regime v — k. 



For a = 0orev = 7r/2 [panels (a) and (b), respectively, 
of Figs. [Hand [5], the two rows of states march past each 
other in an orderly fashion, undergoing a sequence of 
level crossings in complete oblivion of each other's pres- 
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ence. For a = 7r/4 [panel (c) of Figs. [2] and [3], on the 
other hand, all states are part of a coordinated clockwise 
looping motion. While every individual state maintains 
the same position in the level sequence, the wave-shaped 
top row of states has the appearance of moving steadily 
to the right and the bottom row to the left. The path of 
an individual state in the plane of invariants is not un- 
like that of an H2O molecule in a traveling surface water 
wave. 

This qualitative change in pattern caused by differ- 
ent settings of the parameter a requires the assignment 
of mutually exclusive sets of quantum numbers to the 
same set of eigenstates in different parameter regimes. 
The action quantum numbers to, n are the trademark of 
quantum integrability. Their very existence accommo- 
dates level crossings between states of equal parity. The 
level sorting quantum number k, on the other hand, is 
applicable when level crossings between states of equal 
parity are prohibited. It is the trademark of quantum 
nonintegrability. 

This conclusion brings us full circle to the thought 
experiment on invariant tori described in Sec. [U If we 
track an eigenstate along a closed path in (A, a)-plane, 
specifically a path that lies partly inside the integrable 
regime and partly outside, its individuality cannot be 
maintained through a unique and consistent assignment 
of quantum numbers. On a path that first leads a certain 
stretch through the integrable regime and then returns 
through the nonintegrable regime, the tagged eigenstate 
may undergo several crossings on the first leg of this path 
and will then, on the second leg, be unable to cross back 
to its initial position in the level sequence. Barring a mi- 
nor caveat (see Appendix [A"|) this conflict in the assign- 
ment of quantum numbers to eigenstates is a dependable 
detecting device for the demarcation of regimes of inte- 
grability and nonintegrability in quantum systems with 
few degrees of freedom. 



connectedness of the integrable regime in the parameter 
space. Both causes are readily identified as extraneous. 
In the context of the spin-boson model ([T]) only the first 
cause comes into play. 

In the following we describe one scenario where two 
eigenstates swap positions in the level spectrum when 
tracked along a closed path in parameter space, a path 
that does not leave the integrable regime. For this pur- 
pose we consider (TT|) with a = | in the extended param- 
eter space (A,ujs,lub) at a — 0. The energy eigenvalues 



1 



1 



E± = ( n+ - huj B ± -^AA 2 (n + l) + (hw s -hLU B ) 2 , 

(Al) 

and the eigenvectors 

\+) = cos(f>\0,n) +smcf>\l,n), (A2a) 
|-) = -sin0|O,n) +cos0|l,rc), (A2b) 

depend on the angular variable 

E + — nfojJB — 



arctan 



AVnTT 



(A3) 



The point of higher symmetry is at A = 0,Wb = &s- 
Here the energy eigenvalues become doubly degenerate 
(for ujb > 0). We consider the quantum invariant 



(A4) 



(S z )± =±icos20 



defined by expectation values in the eigenstates (|A2[) . 
The loop in parameter space is parametrized as follows: 



hois = Swb(1 + sin 0), 
A = Tiujb (1 — cos 0) , 



(A5a) 
(A5b) 



APPENDIX A: POINT OF HIGHER SYMMETRY 

Conflicts in the assignment of quantum numbers to 
eigenstates may arise for reasons unrelated to noninte- 
grability. In a study of a two-spin system [22| two such 
causes were identified: (i) the presence of points of higher 
symmetry inside the integrable regime; (ii) a multiple 



where < (3 < 2ir. It cuts through the point of higher 
symmetry at = 0. The crucial point is that one com- 
plete loop along this path advances the angle (|A3p by 
Acf> = 7r/2, which interchanges the two states (|A2[) and 
does not bring both invariants (|A1[) and (|A4[) back to the 
same position. It takes two loops to return the states |±) 
to their original identity and the points (E±, (S z )±) to 
their original position. 
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